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The collective electronic excitation in planar sodium clusters is studied by time-dependent density functional 
theory calculations. The formation and development of the resonances in photoabsoiption spectra are inves- 
tigated in terms of the shape and size of the 2-dimensional (2-D) systems. The nature of these resonances 
is revealed by the frequency-resolved induced charge densities present on a real-space grid. For long double 
chains, the excitation is similar to that in long single atomic chains, showing longitudinal modes, end and central 
transverse modes. However, for 2-D planes consisting of (n x n) atoms with n up to 16, new 2-D characteristic 
modes emerge regardless of the symmetries considered. For a kick parallel to the plane, besides the equivalent 
end mode, mixed modes with contrary polarity occur, while for an impulse perpendicular to the plane there will 
be comer, side center, bulk center, and circuit modes. Our calculation reveals the importance of dimensionality 
for plasmon excitation and how it evolves from 1-D to 2-D. 
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INTRODUCTION 

Plasmons are collective electronic excitations in a metal 
caused by an electromagnetic field, which is generally related 
to the surface plasmon resonance (SPR) (lll|2|]. The SPR ef- 
fect in 10-1000 nm wavelength region can be explained as a 
collective oscillation of free electrons responding to the elec- 
tromagnetic perturbation. However, for systems consisting of 
only tens or hundreds of atoms and of a size bridging the 
atomic and condensed phase scales, they feature several ex- 
traordinary properties which are neither exactly the same as 
the Mie plasmon in a classical metal sphere [3] nor like the 
compressional bulk plasmon. The difference is related to the 
existence of discrete energy levels in small nanoparticles or 
atomic clusters, in contrast to the situation of continuous en- 
ergy bands present in solids H. 

The ability of modem technology, like lithography tech- 
niques, to produce nanoparticle arrays with uniform size and 
shape as well as regular inteiparticle separations has made 
it possible to investigate quantitatively the far- and near-field 
interactions of the surface plasmons in these periodic metal- 
nanopaiticle structures 115 - 15 1 On the other hand, some artifi- 
cial atomic structures lll6l - ll9ll with different sizes and shapes 
were also successfully fabricated by making use of the tech- 
nique of scanning tunneling microscope (STM), which have 
been used as ideal prototpye structures for theoretical stud- 
ies of collective electron oscillations (plasmon), because the 
size and shape of these artificial structures are tunable down 
to the precision of a single atom. In particular, the successful 
assembly of linear Au atomic chains on NiAl (110) surface 
with the STM technique jlTIl has arisen an extensive interest 
in studying the plasmon excitations in small 1 -dimensional 
(1-D) atomic chains [4, .20-23.1. 

Theoretically, it is important to investigate the collective 
optical response of these short 1-D systems directly from the 
first principles because Mie's theory pT is not applicable. 



In this regards the time-dependent density functional theory 
(TDDFT) [20.1 has been approved as a powerful tool due to 
its ability of dealing with excited states and its high compu- 
tational efficiency [24]. By using TDDFT calculation in the 
frequency domain Lian et. al. [4] studied the possibility of 
plasmon excitations in linear Au atomic chains with different 
lengths and found that the collective plasmon modes in the ab- 
sorption spectrum arises when the number of atoms is larger 
than ^10. Yuan et al. 121112211 studied the electronic exci- 
tations and the nature of end and central plasmon resonances 
in linear atomic chains of simple and noble metals by using a 
real-space and real-time TDDFT method. They found that in 
case of Ag atomic chains the d electrons can lead to decreased 
intensity and energy of the transverse modes but do not affect 
much the longitudinal modes. Nayyar et al. studied the opti- 
cal properties of short Au atomic chains doped with different 
transition metal atoms (Ni, Rh and Fe) by using the TDDFT 
approach and found that a weak doping with some transition 
metal atoms may lead to the generation of local plasmonic 
modes and, therefore, the transition metal atoms can be used 
to tune the optical properties of nanostructures ll23ll . 



As is well known, dimensionality is one of the most impor- 
tant material-defining parameters: The same chemical com- 
ponents can exhibit dramatically different properties depend- 
ing on whether the atoms are aiTanged in a 0-, 1-, 2- or 3-D 
manner Despite the several recent theoretical studies on the 
plasmon excitations in the 1-D atomic chains, few attention, to 
our best knowledge, has been paid to the plasmonic properties 
of 2-D atomic planes which are now possible to fabricate ex- 
perimentally and may show nontrivial different behavior from 
the coiTesponding 1-D chains. In this work, we investigate 
the electronic excitations of 2-D sodium atomic planes with 
different sizes and shapes by performing TDDFT calculation 
in the time domain. The photoabsoiption spectra due to the 
collective excitations caused by an impulse parallel or per- 
pendicular to the atomic plane are studied. The nature of each 



resonance peak in the spectra is revealed by the frequency- 
resolved induced charge density present on a real-space grid. 
Compared to the related single Na atomic chains which have 
been studied previously as prototype metallic chains l 2Ul22ll . 
the plasmon excitations in these 2-D systems exhibit unique 
characteristic features in the spatial distribution of their in- 
duced charge density, showing some completely different 2-D 
behavior Our result reveals the importance of dimensional- 
ity in the formation and development of plasmon excitations 
in low-dimensional metallic structures and how it will evolve 
from 1-D to 2-D. 



COMPUTATION 

The Na atomic planes studied in the present work are as- 
semblies of the corresponding single Na atomic chains with 
a fixed intra-chain atomic separation and inter-chain spacing 
which are chosen to be a = 2.89 A and b ~ 4.08 A, respec- 
tively, being the same as in the experiment of building linear 
atomic chains on a NiAl surface 11711 where the spacings are 
determined by the lattice constants of the surface. We label 
an atomic plane by {m x n) where m is the number of single 
atomic chains and n is the number of atoms in each chain. The 
ratio of n/m defines the shape of the plane. In this work, we 
consider two cases: double chains (2 x n) with n= 2, 3, 4, 6, 
8, and 18, and square planes (nxn) with n= 2, 3, 4, 5, 6, 7, 8, 
12, 16. All the 2-D structures are placed in the xj;-plane with 
the single atomic chains being along the a;-axis (see the insets 
of Figs[T|and[3li. 

Our calculations are carried out with a real-space and real- 
time TDDFT program Octopus\Z5], in which the Na ion is 
described by norm-conserving pseudopotentials and the lo- 
cal density approxiamtion (LDA) is used for the electron ex- 
change and correlation for both the ground-state and excited- 
state calculations. A grid in real space which is defined by 
assigning a sphere around each atom with a radius of 6 A 




FIG. 1 : The dipole response of sodium double atomic chains (2 x 
n) with a varying length n. The interatomic distances in the x and 
y directions are a=2.89A and b=4.08A, respectively. The series of 
spectra correspond, counting from the bottom, to n = 2, 3, 4, 6, 8, 
and 18, respectively. 



and a uniform mesh grid of 0.3 A is adopted to describe the 
wavefunction and charge density. To obtain the excitation 
spectrum, the system is impulsed from its initial ground state 
with a very short delta-function-like perturbation, and then the 
time-dependent Kohn-Sham equation is evolved in real space 
and real time for a certain period of time. Specifically, an 
electronic wave packet is evolved for 10,000 time steps with 
each being 0.003 h/eW long. After the real-time propagation, 
the photoabsorption spectrum is extracted by Fourier trans- 
forming the time-dependent dipole strength. Furthermore, a 
3-D image of the frequency-resolved induced density distri- 
bution in real space is then obtained for each resonance in the 
spectrum by Fourier transforming the time series of the total 
induced charge density at each resonance frequency for ev- 
ery real-space mesh grid and then is visualized by a graphic 

lii. 
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RESULT AND DISCUSSION 



Double chains 



Let us start with the double chains. Fig.[T]shows the evolu- 
tion of the (optical absorption) dipole strength as a function of 
the length (i.e., the number of atoms in each single chain (n)). 
One can see that as the double chain becomes long enough 
(n > 6), its two transverse modes gradually split into four 
branches with two low-energy resonances impulsed along y- 
axis and two high-energy resonances caused by a kick along z- 
axis (perpendicular to the double-chain plane). That is to say, 
compared with the results of single Na chains given by a previ- 
ous calculation (Fig. 1 of Ref. 1.21.1 ). one transverse mode of a 
single chain now split into two due to the asymmetry between 
the y and z directions, but with similar evolution in energy as 
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FIG. 2: Fourier transforms of the induced densities at the frequencies 
of the plasmon resonances indicated on the left for the double chain 
(2 X 18), as shown in Fig.[T] Green and red isosurfaces denote the 
negative and positive densities, respectively. 
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FIG. 3: The dipole response of the asymmetrical square sodium 
planes (n x n) with a varying size n. The interatomic distances in 
the X and y directions are a=2.89A and 6=4. OSA, respectively. The 
series of spectra correspond, counting from the bottom, to n =2, 3, 
4,5, 6,7, 8,12, and 16. 

the length of the chain increases. The low-lying longitudinal 
Lxi mode behaviors similarly as in the case of single chain: It 
redshifts in energy as the double chain becomes longer, which 
results from the accumulation of collectivity in the dipole os- 
cillation. Its intensity increases near linearly with the chain 
length, which can arise from the reduction of the energy gaps 
involved in the dipole excitation ll22ll . 

The nature of these resonances can be understood by look- 
ing at the induced charge densities at the resonance frequen- 
cies. Fig. |2] shows the frequency-resolved induced densities 
for the double chain with n=18. It can be seen that Lxi is a 
single longitudinal mode and Lx2 is a multipolar longitudinal 
mode which may be resulted from the longitudinal quantiza- 
tion of plasmon, as previously found in long single chains 112211 
. On the other hand, the T^^ (TEy) and TC^ (TCy) peaks 
in Fig. [U correspond to the end and central transverse modes, 
respectively, when the kick is along z-axis (y-axis). The fre- 
quencies of TEy and TCy modes are lower than those of the 
TEz and TCz modes, respectively, because one more atom 
is added in the y direction. In both cases, the TE mode has a 
lower excitation energy than the TC mode. Such a spatial sep- 
aration of the transverse excitations resembles the surface and 
bulkplasmons of solids and thin films, where the surface (end) 
plasmons have lower energies than the bulk (central) plasmon. 
Overall, our calculation shows that the collective excitation in 
long double chains (with a large length-width ratio) is similar 
to that in long single chains although the transverse modes are 
now split. 



Asymmetrical square planes 

Next we consider the case of square planes, (nxn), with dif- 
ferent sizes n = 2, 3, 4, 5, 6, 7, 8, 12, 16, and plot their excita- 
tion spectra in Fig.|3] It is obvious that starting from the same 
smallest cluster (2x2) the evolution of the resonance peaks in 
energy for the square planes is totally different from that for 
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FIG. 4: The excitation energy (a) and dipole strength (b) as a function 
of the size n of the asymmetrical square sodium planes (n x n), for 
the resonance peaks shown in Fig.|3] 



the single or double chains. When the impulse is along the 
X- or y-axis the resulting low-energy resonances Xi and Yi 
are split because of the asymmetry (a y^ 6). One should note 
that the increase in the interatomic separation (b > a) will 
lower the coiTesponding resonance frequency (li is lower in 
energy than Xi). This can be understood by considering the 
fact that the nonlocal s electrons can travel a bit longer dis- 
tance with a larger interatomic separation in the y direction. 
When the square plane becomes large enough (n > 8) some 
other higher-energy resonances (X2, Y2, X3, and Y3 in Fig.O 
begin to emerge. These resonances behave differently from 
Xi and Yi : Now the Y peaks are higher in energy than the 
X peaks, indicating their oscillation modes are different (see 
discussion later). A common trend for all these plasmon res- 
onances in the larger planes is that their energy redshifts and 
their strength increases near linearly with the increased size 
{n) of the plane (see Fig. |4]l. The linear redshift can be un- 
derstood by considering the reduction of the energy gap as the 
cluster becomes larger The increase in strength can be as- 
cribed to the accumulation of collectivity in the excitation as 
more electrons participate in the collective oscillation. 

When the dipole excitation is along the z-axis, there is only 
one mode (Zi) visible for the two smallest clusters (i.e., (2x2) 
and (3 x 3), see Fig. |3]l. As the atomic plane becomes larger, 
more high-energy modes gradually emerge. The first one is 
Z2 and it dominates in intensity starting from (6 x 6). When 
the plane is larger than (8x8) this resonance splits into two 




FIG. 5: Fourier transforms of the induced densities at the frequencies 
of the plasmon resonances excited along the x and y directions for 
the asymmetrical square plane (16 x 16), as shown in Fig.[3l 



peaks (Z2-1 and Z2-2) and then the Z2-2 increases its in- 
tensity rapidly and develops into a dominant plasmon reso- 
nance (a 2-D central mode, as discussed later). Starting from 
(5x5) another high-energy mode (Z3) emerges. As shown 
in Fig. |4] the intensities of all these Z plasmon modes in the 
larger planes increase near linearly as the size of the plane gets 
larger because of the accumulation of collectivity. However, 
the change in their excitation energies is much smaller com- 
pared to the case of the X and Y modes because the dipole 
oscillation is perpendicular to the atomic plane and therefore 
is affected much less by the size of the plane, as similarly 
found in the case of single atomic chains 112711 . 

To understand the nature of the X and Y plasmon res- 
onances in the square Na planes we plot their frequency- 
resolved induced charge densities in Fig. |5]for the (16 x 16) 
system. For the Xi and Yi resonances, the induced density 
is mainly distributed around the four corners of the plane, 
which is somewhat similar to the TE mode in single chains. 
This is understandable because the 2-D plane can be viewed 
as a zonal chain along the direction perpendicular to the im- 
pulse direction. However, the asymmetry between the x- and 
y-direction leads to slightly different excitation energies be- 
tween the Xi and Yi modes as shown in Fig. |3] The Yi has 
a slightly lower energy because the nonlocal s electrons can 
travel a bit longer distance in the y direction. Differing from 
the Xi and Yi modes, the X2, Y2, X3, and Y3 modes have 
a mixed density distribution. For X2, besides the equivalent 
center (EC) mode, an equivalent end (EE) mode with the op- 
posite polarity also exists simultaneously. For Y2, there is a 
EC mode and a multipolar end (ME) mode with the opposite 
polarity existing simultaneously. Similarly, X^ is EEh-ME and 
Y3 is MEh-ME. Because the two simultaneous modes have op- 
posite polarity the total dipole strength of these mixed modes 
is much smaller than Xi and Yi . The different plasmon ex- 
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FIG. 6: Fourier transforms of the induced densities at the frequencies 
of the plasmon resonances excited along the z direction for (a) (5x5), 
(b) (6x6), and (c) (16x 16) asymmetrical square planes, as shown in 

Fig.m 

citation also lead to a different order of their excitation ener- 
gies, for instance, Y2 has a higher energy than X2- Note that 
the asymmetry between the x- and y-direction is the reason 
for the nonequivalence between the X and Y plasmon reso- 
nances, as can be seen clearly in Fig. |5] upper panel vs lower 
panel. 

When an impulse perpendicular to the Na plane is imposed, 
the frequency-resolved induced charge densities of the result- 
ing Z plasmon modes are plotted in Fig. |6] for three sizes of 
the plane ((5 x 5), (6 x 6), and (16 x 16)) to show the evolution 
of these modes. One can see that when the Na plane is small 
((5x5) and (6x6)) the density distribution is quite iiTegularbut 
evolves to form regular patterns when the Na plane becomes 
large enough (16 x 16): Zi turns out to be an equivalent end 
mode with its density located around the longer sides of the 
plane; Z2 evolves gradually to form two modes with their den- 
sity distributed around the center of the shorter sides {Z2-1) 
and in the bulk center (.^2-2); ^3 evolves eventually into a 
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FIG. 7: The dipole response of the symmetrical square planes (n x 
n) with a varying size n. The interatomic distances in the x and 
y directions are a = 6=2. 89 A, respectively. The series of spectra 
correspond, counting from the bottom, to n =2, 3, 4,5, 6,7, 8 ,12, 
and 16. 
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FIG. 8: Fourier transforms of the induced densities at the frequencies 
of the plasmon resonances excited along the x and y directions for 
the symmetrical (16 x 16) square plane, as shown in Fig. [7] 

circuit mode. Like in the case of single and double atomic 
chains, here the end/surface modes (Zi and Z2-1) have lower 
excitation energies than the central/bulk mode (^2-2)- Note 
that the induced density distributions from all the Z plasmon 
modes are highly complementary in space. 
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FIG. 9: Fourier transforms of the induced densities at the frequencies 
of the plasmon resonances excited along the z direction for (a) (5x5), 
(b) (6 X 6), and (c) (16 x 16) symmetrical square planes, as shown in 
Fig.H 



Symmetrical square planes 

From the density distribution of the Zi or Z2-1 plasmon 
mode one can see clearly the significant effect from the asym- 
metry of the plane (a ^ b). To have a clearer picture 
of this effect and to see more generally the plasmonic be- 
havior in 2-D systems, we also perform calculations for the 
symmetrical Na planes with a = b = 2.89 A. The photoab- 
sorption spectra are given in Fig. [Tjfor the systems (n x n), 
ri = 2, 3, 4, 5, 6, 7, 8, 12, 16. By comparing these results with 
those of the asymmetrical planes (Fig. |3] ) one can see that 
now the X and corresponding Y peaks merge into single ones 
(i.e., Xi, X2, and X^ peaks in Fig.|7]i because of the symme- 
try. The energy of the Xi (symmetrical) is also quite close to 
the Xi (asymmetrical) because now lattice constants in both 
X and y directions are the same as the lattice constant in the 
X direction of the asymmetrical systems. By looking at the 
frequency -resolved induced densities shown in Fig.|8]one can 
also see that the density distributions of the Xi 2.3 modes in 



the symmetrical (16 x 16) plane are very close to those of the 
-^^1,2.3 in the asymmetrical one but are remarkably different 
from those of the Yi.2,3. 

On the other hand, this change in symmetry does not, how- 
ever, affect qualitatively the Z modes: the overall shape of 
the spectra almost remains the same (Fig. Q vs Fig. [Sj ex- 
cept for the remarkable blueshifts in Fig. Qand a new high- 
energy resonance (^3-2) emerging in (16 x 16). To see the 
nature of these remarkably blueshifted Z modes, we plot their 
frequency-resolved induced density in Fig. |9]for the symmet- 
rical (5x5), (6x6), (16x16) planes. One can see that now the 
charge distribution is symmetrical in the x and y directions 
for all the Z resonances because of the symmetry which also 
leads to a quite different evolution and the eventual pattern in 
the large(16xl6)plane for each Zmode. Now Zi evolves into 
a pure corner mode instead of the longer-side mode. Z2 de- 
velops finally into two modes, a side-center mode and a strong 
bulk center mode which is similar to that in the asymmetrical 
case and is the most dominant mode in dipole strength in both 
cases. Z3 evolves eventually into a circuit mode and a weak 



central mode. Like in the case of asymmetry, here the induced 
densities from the different modes are highly complementary 
in space. 

Our calculations show that the plasmon excitation in 2-D 
atomic structures is sensitive to the structural symmetry. In 
the present case, a change of symmetry in the xj/-plane will 
affect qualitatively the X and Y resonances but will not have 
a qualitative effect on the Z resonances although the quanti- 
tative effects can be significant, like the remarkable blueshifts 
and the changes in the detailed pattern of collective charge 
oscillation. 



SUMMARY 

We have studied the collective electronic excitations in pla- 
nar sodium clusters by performing the time-dependent density 
functional theory calculations in the time domain. Two kinds 
of systems are considered: double chains and square atomic 
planes. The formation and development of the collective res- 
onances in their photoabsorption spectra are investigated as a 
function of their size. The nature of these plasmon modes is 
revealed by the frequency-resolved induced charge densities 
present on a real-space grid. 

For the double chains with a large length-width ratio, the 
longitudinal modes are similar to those in single chains, while 
the transverse modes show two pairs of bimodal structure be- 
cause of the lower symmetry, which can be still assigned as the 
end and central modes, respectively, just as in the case of sin- 
gle atomic chains. However, for larger square atomic planes 
consisting of {n x n) atoms with n up to 16, new 2-D char- 
acteristic modes emerge. For a kick parallel to the plane, be- 
sides the equivalent end mode, the induced density also shows 
mixed modes with contrary polarity. On the other hand, for an 
impulse perpendicular to the plane the induced density shows 
2-D characteristic corner, side center, bulk center, and circuit 
modes. 

The structural symmetry of an atomic plane plays a signif- 
icant role in its plasmon excitation. Changing a atomic plane 
from being in-plane symmetrical to asymmetrical will lead to 
splitting of the degenerated in-plane modes. It will also cause 
remarkable changes in the nature and excitation energies of 
the perpendicular modes, but will not change qualitatively the 
overall shape of the perpendicular photoabsorption spectra. 

Our calculations reveal the importance of dimensionality in 
the formation and development of plasmon excitation in low- 
dimensional metallic structures and how it will evolve from 
l-Dto2-D. 
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